What is claimed is: 



1 . A method for blind source separation of multiple sources, comprising: 
detecting the multiple sources using an array of sensors to obtain data 

representative of the multiple sources; 

representing the data by two mixtures having estimates of amplitude and delay 

mixing parameters; 

updating the estimates of amplitude and delay mixing parameters, comprising the 
steps of: 

calculating a plurality of error measures, each of the plurality of error 
measures indicating a closeness of the estimates of amplitude and delay mixing 
parameters for a given source to a given time-frequency point in the two mixtures; and 

revising the estimates of amplitude and delay mixing parameters, based on 
the plurality of error measures; 

filtering the two mixtures to obtain estimates of the multiple sources, comprising 
the steps of: 

selecting one of the plurality of error measures having a smallest value in 
relation to any other of the plurality of error measures, for each of a plurality of time- 
frequency points in the mixtures; and 

leaving unaltered any of the plurality of time-frequency points in the 
mixtures for which a given one of the plurality of error measures has the smallest value, 
while setting to zero any other of the plurality of time-frequency points in the mixtures 
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for which the given one of the plurality of error measures does not have the smallest 
value, for each of the plurality of error measures; 

outputting the estimates of the multiple sources. 

2. The method of claim 1 , wherein said step of representing the data by two 
mixtures comprises the step of computing a first mixture xl and a second mixture x2 as 
follows: 
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where N is a number of the multiple sources; 6j is an arrival delay between the 
array of sensors resulting from an angle of arrival; aj is a relative attenuation factor 
corresponding to a ratio of attenuations of paths between the multiple sources and the 
array of sensors; Sj (t) is a j th source; j is a source index ranging from 1 to N, where N is a 

number of the multiple sources; and t is a time argument. 

3 . The method of claim 1 , wherein said step of calculating the plurality of 
error measures comprises the step of computing 
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where p(a j9 S j9 a) 9 T k ) is an error measure for a j l source; aj and 8j are current 
estimates of amplitude and delay mixing parameters, respectively; X x (co,T k ) and 
X 2 (co 3 T k ) are time-frequency representations of a first mixture and a second mixture of 
the two mixtures, respectively; k is a current time index; x k is a time argument 
corresponding to a k th time index; co is a frequency argument; and j is a source index 
ranging from 1 to N, where N is a number of the multiple sources. 

4. The method of claim 1, wherein said step of revising the estimates of 
amplitude and delay mixing parameters comprises the step of computing a magnitude and 
a direction of change for each of the estimates of amplitude and delay mixing parameters. 

5. The method of claim 4, wherein said step of computing the direction of 
change for each of the estimates of amplitude mixing parameters comprises the step of 
computing 
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where ^ J ^ x ^ represents the magnitude and the direction of change in a current 

da } 

estimate of amplitude mixing parameter for a j* source that causes a largest change in 
correspondence between the current estimate and the data; k is a current time index; x k is 
a time argument corresponding to a k th time index; aj is the current estimate of amplitude 
mixing parameter for the j th source; Sj is a current estimate of delay mixing parameter for 
the j th source; X x (a>,T k ) and X 2 (co ? x k ) are time-frequency representations of a first 
mixture and a second mixture of the two mixtures, respectively; o is a frequency 
argument; j and / are source indexes ranging from 1 to N, where N is a number of the 
multiple sources; p(a J9 5 j9 a) 9 T k ) is an error measure for the j th source; A, is a smoothness 
parameter; and Re is a function that returns a real part of a complex number. 

6. The method of claim 4, wherein said step of computing the direction of 
change for each of the estimates of delay mixing parameters comprises the step of 
computing 
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where represents the magnitude and the direction of change in a current 

68 j 

estimate of delay mixing parameter for a j th source that causes a largest change in 
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correspondence between the current estimate and the data; k is a current time index; x k is 
a time argument corresponding to a k th time index; aj is a current estimate of amplitude 
mixing parameter for the j th source; 85 is the current estimate of delay mixing parameter 
for the j th source; X x (G>,T k ) and X 2 (a>,x k ) are time-frequency representations of a first 
mixture and a second mixture of the two mixtures, respectively; © is a frequency 
argument; j and / are source indexes ranging from 1 to N, where N is a number of the 
multiple sources; p{a j ,8 j ,a> y t k ) is an error measure for the j th source; A, is a smoothness 

parameter; and Im is a function that returns an imaginary part of a complex number. 

7. The method of claim 4, wherein said step of revising the estimates of 
amplitude and delay mixing parameters comprises the step of computing 



aj [k] = aj[k-l]-fiaj[k]^^ 



8j[k) = Sjik-Vl-pSjlk]-^ 



where aj[k] and 85 [k] are the estimates of amplitude and delay mixing parameters 

for a j source at a time index k, respectively; B is a learning rate constant; — 

da } 

represents the magnitude and the direction of change in a current estimate of amplitude 
mixing parameter for a j th source that causes a largest change in correspondence between 
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the current estimate and the data; — — represents the magnitude and the direction of 

change in a current estimate of delay mixing parameter for a j th source that causes a 
largest change in correspondence between the current estimate and the data; r k is a 
current time; k is a current time index; and j is a source index ranging from 1 to N, 
wherein N is a number of the multiple sources. 

8. The method of claim 4, wherein said step of computing the magnitude of 
change for each of the estimates of amplitude and delay mixing parameters comprises the 
step of scaling the magnitude of change by a learning rate constant and a variable learning 
rate. 

9. The method of claim 8, wherein said step of computing the learning rate 
constant comprises the step of computing 

I 

where qy[k] represents an amount of mixture energy that is defined by estimates of 
amplitude and delay mixing parameters for a j th source; k is a current time index; x k is a 
time argument corresponding to a k th time index; a s and 8j are current estimates of 
amplitude and delay mixing parameters, respectively for the j th source; ^(co/Ck) and 
X 2 (co^ik) are time-frequency representations of a first mixture and a second mixture of 
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the two mixtures, respectively; co is a frequency argument; j and / are source indexes 
ranging from 1 to N, where N is a number of the multiple sources; and X is a smoothness 
parameter. 



10. The method of claim 8, wherein said step of computing the variable 
learning rate comprises the step of computing 



mixture energy that is defined by estimates of amplitude and delay mixing parameters for 
a j tb source; y is a forgetting factor; m is a time index ranging from 0 to a current time 
index k; and j is a source index ranging from 1 to N, where N is a number of the multiple 
sources. 

1 1 . The method of claim 1 , wherein said step of selecting one of the plurality 
of measures comprises the step of computing 



CCj[k] = 



qj[k] 



where a } [k] represents the variable learning rate; qj[k] represents an amount of 





0 otherwise 
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where Qj(co,x k ) is a time-frequency mask; p(a J9 S j9 a> 9 T k ) is an error measure for a j 

source; j and m are source indexes ranging from 1 to N, where N is a number of the 
multiple sources; v k is a current time; Hs a current time index; aj and 5j are current 
estimates of amplitude and delay mixing parameters for the j th source, respectively; and co 
is a frequency argument. 

12. The method of claim 1 , wherein said leaving and setting steps comprise 
the step of computing 

Sj (Q) 9 T k ) = Q J (0),T k )X l (0> 9 T k ) 

where Sj(co,x k ) is an estimate of a time-frequency representation of a j th source; j 
is a source index ranging from 1 to N, where N is a number of the multiple sources; 
Qj(co ? x k ) is a time-frequency mask; X x (oo ? T k ) is a time-frequency representation of a first 
mixture of the two mixtures; co is a frequency argument; r k is a current time; and A: is a 
current time index. 

1 3 . The method of claim 1 , wherein said step of selecting one of the plurality 
of measures comprise the step of constructing a time-frequency mask based on a 
maximum likelihood parameter estimation. 
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14. The method of claim 13, wherein said step of constructing the time- 
frequency mask comprises the step of computing 



measure for a j th source; j and m are source indexes ranging from 1 to N, wherein N is a 
number of the multiple sources; r k is a current time; k is a current time index; a } and 5j are 
current estimates of amplitude and delay mixing parameters for the j th source, 
respectively; and co is a frequency argument. 

15. The method of claim 1 , wherein said leaving and setting steps comprise 
the step of computing estimates for time-frequency representations of each of the multiple 
sources. 

1 6. The method of claim 1 5, wherein said step of computing the estimates for 
time-frequency representations comprises the step of computing an estimate for a j th 
source of the multiple sources as follows: 




where Qj(co,Tk) is a time-frequency mask; p{a } , 8 j , a>, r k ) is an error 



Sj (a> 9 r k ) = Q,. (o>, t k ) X x r k ) 



where Sj(co,T k ) is an estimate of a time-frequency representation of a j source; j 
is a source index ranging from 1 to N, where N is a number of the multiple sources; 
Qj(co ? Tk) is a time-frequency mask; X x (co,T k ) is a time-frequency representation of a first 
mixture of the two mixtures; o is a frequency argument; r k is a current time; k is a current 
time index. 

17. The method of claim 1, wherein said filtering step comprises the step of 
partitioning a time-frequency representation of one of the two mixtures to demix all of the 
multiple sources. 

1 8. The method of claim 1, wherein said outputting step comprises the step of 
reconstructing the multiple sources using a dual window function applied to the 
estimates. 

1 9. The method of claim 1 , further comprising the step of computing a 
frequency domain representation for the two mixtures, prior to said updating step. 

20. The method of claim 1 , further comprising the step of computing a time 
domain representation for the estimates of the multiple sources, prior to said outputting 
step. 
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21 . The method of claim 1 , wherein said method is implemented by a program 
storage device readable by machine, tangibly embodying a program of instructions 
executable by the machine to perform said method steps. 
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